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by 
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SUMMARY 

Chemical  nucleation  has  been  studied  numerically  in  a  stagnation  point  mixing  layer 
in  which  reactants  in  two  counter-flowing  streams  form  a  condensable  monomer.  The 
response  of  the  subseq^.ent  nucleation  kinetics  to  the  velocity  gradient  in  the  flow  is 
described  in  terms  of  a  Damkohler  number.  Two  limiting  cases  have  been  established. 
Firstly,  if  the  Damkohler  number  for  monomer  production  is  .small,  i.e.  the  rate  of  monomer 
production  is  slow,  then  the  nucleation  of  particles  can  be  strongly  affected  by  the  flow 
field  in  a  manner  which  is  equivalent  to  the  effect  of  supersaturation  in  a  uniform  vapour. 
Secondly,  if  the  Damkohler  numbers  for  clu.ster  growth  are  small  (due  to  a  small  accommo¬ 
dation  factor  for  monomer-cluster  interactions),  the  concentrations  of  clusters  do  not 
achieve  equilibrium  levels.  This  can  re.sult  in  the  suppression  of  particle  formation  over  a 
critical  range  of  Damkohler  numbers.  In  this  case  the  behaviour  of  the  nucleation  kinetics 
is  analogous  to  the  tran.sient  phase  of  nucleation  in  a  uniform  vapour. 
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NOMENCLATURE 


velocity  gradient 
cluster 

multicomponent  diifusivity 
Damkohler  number 
internal  energy 
modified  stream  function 
Gibbs  free  energy 
standardised  enthalpy 
nucleation  rate 

Boltzmann’s  constant  (equations  2-4) 
rate  constant 

number  of  cluster  size  groups 
molecular  weight 
number  of  monomers 
pressure 

gas  constant  [m®  atm  mole“i  K“i] 

gas  constant 

radius 

standardised  entropy 

supersaturation  ratio 

Schmidt  number 

temperature 

tangential  velocity 

volume  of  molecule  in  liquid  phase 

normal  velocity 

molar  rate  of  production 

tangential  coordinate 

normal  coordinate 

mass  fraction 

mass  flux 

mass  accommodation  factor 
non-dimensionalised  normal  coordinate 


Superscripts 


dynamic  viscosity 
kinematic  viscosity 
density 

surface  tension 


critical  cluster 


differentiation  with  respect  to  ^ 
differentiation  with  respect  to  time 
at  standard  conditions 


Subscripts 


reactant  A 
reactant  B 

at  edge  of  mixing  layer 
at  thermodynamic  equilibrium  (Section  2.1) 
inert  gas 

refers  to  species  i  with  i  monomer  units 
refers  to  species  j  with  j  monomer  units 
liquid  phase 

number  of  monomer  units 
product 


CAPTIONS  TO  nCURES 


FIGURE  I  Variation  in  free  energy  functions  with  cluster  size. 

FIGURE  2  Stagnation  point  mixing  layer  flow  used  in  the  calculations. 

FIGURE  3  Mixing  layer  velocity  profiles. 

FIGURE  4  Mass  fractions  of  reactants  A  and  B  with  distance  from  cylinder. 

FIGURE  5  Estimated  Gibbs  free  energy  functions  for  monomer  association. 

FIGURE  6  Cluster  mass  fractions  at  stagnation  point  C  =  10“'*  and  AGn  for  NH4CI  at  400  K. 

FIGURE  7  Mass  fractions  and  concentrations  at  stagnation  point  of  particles  containing 

200  monomer  units 

(  =  10-*  and  AGn  for  NH4CI  at  400  K. 

FIGURE  8  Distribution  of  clusters  through  the  flow  for  clusters  with 

(a)  2  monomer  units 

(b)  40  monomer  units 

(c)  80  monomer  units 

?  =  10-4  and  AGn  for  NH4CI  at  400  K. 

FIGURE  9  Cluster  mass  fractions  at  stagnation  point 

{  =  10-4  and  AGn  for  NH4CI  at  550  K. 

FIGURE  10  Cluster  mass  fractions  at  stagnation  point  for  slow  monomer  formation  rates 
S  =  1  and  AGn  for  NH4CI  at  400  K. 

FIGURE  11  Mass  fractions  and  concentrations  at  stagnation  point  of  clusters  containing 
200  monomer  units 

^  =  1  and  AGn  for  NH4CI  at  400  K. 

FIGURE  12  Cluster  mass  fractions  at  stagnation  point  for  slow  monomer  formation  rates 
C  =  1  and  AGn  for  MgO  at  1700  K. 


1.  INTRODUCTION 


Condensation  of  a  vapour  to  a  liquid  in  the  absence  of  foreign  particles  can  be  described  by 
homogeneous  nucleation  theory.  Experiments  usually  involve  raising  the  supersaturation  of  the 
condensable  vapour  to  a  critical  level  at  which  observable  liquid  or  solid  phase  particles  are 
formed.  The  nucleation  rate  which  is  determined  by  such  experiments  refers  to  the  formation 
rate  of  critically  sized  clusters  of  vapour  molecules  (known  as  nuclei)  which  tend  to  grow  as 
quickly  as  they  evaporate. 

In  the  classical  theory  the  formation  of  a  critical  cluster  is  achieved  by  the  addition  of 
monomer  or  single  vapour  molecules  to  sub-critical  clusters.  The  tendency  of  a  cluster  to  re¬ 
evaporate  is  determined  by  a  free  energy  function.  Evaluation  of  this  free  energy  function  can 
be  achieved  by  use  of  a  so-called  liquid  droplet  model  in  which  the  properties  of  a  growing 
cluster  are  identified  with  the  properties  of  the  bulk  liquid  phase.  Later  developments  have 
sought  to  improve  on  this  estimation  of  the  free  energy  function.  Springer  (1978)  has  reviewed 
recent  advances  in  homogeneous  nucleation  theory. 

Homogeneous  nucleation  does  not  always  proceed  simply  by  the  supersaturation  of  pure, 
condensable  vapour.  Chemical  nucleation  designates  a  process  in  which  the  initial  formation 
of  a  condensable  vapour  or  its  subsequent  nucleation  proceeds  by  means  of  chemical  reactions. 
Formation  of  a  condensable  vapour  by  chemical  reaction  can  occur  in  combustion  systems 
such  as  rocket  exhausts  (Courtney,  1962),  or  in  the  formation  of  inorganic  oxide  aerosols  during 
coal  combustion  (Katz  and  Donohue  (1982),  Neville  and  Sarofim  (1982)).  The  production  of 
atmospheric  aerosols  may  also  be  the  result  of  chemical  reaction  (Friedlander,  1977).  The 
nucleation  of  soot  particles  in  a  flame  can  be  described  as  a  chemical  nucleation  process  in  the 
sense  that  small  molecules  add  to  molecules  with  increasingly  large  carbon  numbers  in  a  chain 
of  reactions  that  leads  to  a  macroscopic  solid  particle.  In  contrast  to  physical  condensation 
each  step  involves  chemical  bonding:  clusters  arc  no  longer  /t-mers  formed  from  a  single 
monomer  species  but  are  distinct  chemical  compounds. 

Most  practical  situations  in  which  chemical  nucleation  occurs  will  involve  mixing  of  vapour- 
phase  reactants.  Mixing  will  proceed  by  dilfusive  and  convective  transport. 

Diffusive  convection  in  chemically  reacting  flows  has  been  analysed  by  Fendell  (1965)  who 
considered  the  effect  of  the  straining  motion  in  a  stagnation  point  flow  on  ignition  and  extinction 
of  a  diffusion  flame.  The  interaction  between  diffusion,  convection  and  reaction  was  analysed 
in  terms  of  a  Damkohler  number  which  was  the  ratio  between  a  characteristic  flow  time  scale  and 
a  characteristic  chemical  reaction  lime  scale.  Extinction  of  a  flame  was  predicted  to  occur  at  a 
sufficiently  low  critical  Damkohler  number,  i.e.  when  the  velocity  gradient  in  the  stagnation 
flow  was  sufficiently  high.  Values  of  the  critical  velocity  gradient  for  counter-flow  diffusion 
flames  with  various  fuels  have  been  measured  (Tsuji,  1982).  These  theoretical  and  experimental 
studies  were  concerned  with  high  activation  energy  flame  reactions  which  exhibited  extinction 
phenomena  in  a  stagnation  point  flow.  They  do  suggest,  however,  a  role  for  the  flow  yield  in 
influencing  the  essentially  isothermal  reactions  involved  in  chemical  nucleation  in  a  stagnation 
point  mixing  layer. 

The  involvement  of  transport  processes  in  nucleation  introduces  additional  phenomena 
which  are  not  considered  in  the  classical  theory  for  a  uniform  vapour  (for  example,  Becker  and 
Reiss,  1976).  For  most  simple  substances,  such  as  metal  vapour  which  do  not  involve  chemical 
reactions  during  their  condensation,  the  time  scale  for  cluster  growth  and  nucleation  is  con¬ 
siderably  shorter  than  the  residence  time  of  the  condensing  vapour  in  a  flow.  This  situation 
applies  even  in  high  speed  nozzle  flows  (Wegener  and  Wu,  1977).  However,  if  the  time  scale  for 
nucleation  processes  is  increased  significantly  by  chemical  reaction  so  that  it  is  comparable  to  a 
characteristic  time  scale  for  the  mixing  flow,  the  flow  may  affect  the  finite  rate  nucleation  kinetics. 
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The  classical  theory  of  nucleation  is  concerned  with  a  steady-state  process.  The  transient 
phenomenon  which  occurs  between  the  achievement  of  a  critical  supersaturation  and  the  onset 
of  a  steady  nucleation  rate  has  been  studied  numerically  in  a  condensing,  uniformly  distributed 
water  vapour  by  Courtney  (1962)  and  Abraham  (1969).  Analytical  treatments  of  nucleation 
kinetics  by  Wakeshima  (1954),  Collins  (1955),  Andres  and  Boudart  (1965)  and  others  are  essen¬ 
tially  in  agreement  with  the  numerical  predictions  for  water  condensation.  The  time  lag  to  a 
steady  state  is  of  the  order  of  1  /i  sec;  this  is  considerably  shorter  than  the  characteristic 
time  scale  for  most  practical  flows.  The  analysis  of  the  time  lag  shows,  however,  a  direct  pro¬ 
portionality  to  a  mass  accommodation  or  sticking  factor  which  is  an  indication  of  the  fraction 
of  monomer-cluster  collisions  that  are  successful  in  increasing  the  cluster  size.  This  factor  is 
believed  to  be  of  order  unity  for  substances  such  as  water;  for  substances  which  involve 
chemical  reaction  in  their  nucleation  it  may  be  much  less  than  one.  In  the  context  of  nucleation 
kinetics,  an  activation  energy  (£)  for  the  reaction  leading  to  condensable  monomer  or  for  the 
subsequent  monomer-cluster  reactions  may  be  associated  with  an  accommodation  factor, 
i.e.  the  term  e  e  rt  an  Arrehnius  rate  expression  is  equivalent  to  the  accommodation  or 
sticking  factor  commonly  applied  in  physical  condensation  processes. 

The  present  study  was  prompted  by  the  observation  of  soot  or  carbon  particle  formation 
in  a  gaseous  diffusion  flame.  The  presence  of  the  soot  aerosol  accounts  for  the  yellow  luminosity 
observed  in  a  sooty  flame.  Tsuji  and  Yamaoka  (1968.  1970)  have  reported  gas  sampling  measure¬ 
ments  in  laminar  counter-flow  diffusion  flames  that  were  stabilized  in  the  forward  stagnation 
point  of  a  porous  cylinder.  Fuel  was  blown  from  the  cylinder  into  a  stream  of  air.  They  reported 
a  critical  velocity  gradient  (approach  flow  velocity  divided  by  burner  radius)  above  which  no 
soot  was  observed  to  form  in  the  flame,  regardless  of  the  rate  at  which  fuel  was  blown  from  the 
burner.  For  propane  and  air  this  critical  velocity  gradient  was  about  100  s  *. 

Disappearance  of  the  soot  occurs  at  a  velocity  gradient  which  is  substantially  less  that  the 
velocity  gradient  associated  with  extinction  of  the  flame  itself  (Tsuji,  1982).  Flame  temperature 
is  quite  insensitive  to  the  approach  flow  over  the  velocity  range  in  which  soot  disappears  and, 
consequently,  it  is  unlikely  that  the  effect  has  a  chemical  origin.  The  sensitivity  of  the  soot 
formation  mechanism  to  the  approach  flow  velocity  is  quite  marked  in  this  flame.  A  velocity 
gradient  somewhat  less  than  the  critical  condition  yields  a  large  number  concentration  of  particles 
in  the  flame,  typically  of  the  order  of  10’®  to  10’"  m  ”  (Vandsburger,  Kennedy  and  Glassman, 
(1984)).  An  increase  in  the  velocity  gradient  through  increasing  the  approach  flow'  velocity  by 
as  little  as  50%  can  reduce  the  particle  number  concentration  to  zero. 

The  formation  or  nucleation  of  the  initial,  solid  soot  particles  involves  a  complex  chemical 
process  which  is  the  subject  of  continuing  debate  (Cilassman,  Brezinsky.  Gomez  and 
Takahaski  (1984),  Calcotc  (1981),  Lahayc  and  Prado  (1974)).  Regardless  of  the  details  of  the 
particular  mechanism  it  is  apparent  that  the  formation  of  a  soot  particle  involves  the  addition 
of  small  molecules  to  growing  larger  molecules  in  a  scries  ol' reactions  leading  up  to  a  macroscopic 
particle.  The  addition  of  each  small  molecule  involves  chemical  bonding  ofteii  with  an  associated 
activation  energy  which  results  in  a  slow  growth  rale  compared  with  the  collision  rate.  It  is  one 
of  the  aims  of  this  study  to  explore  the  possible  interaction  between  the  flow  field  and  the  kinetics 
of  a  nucleation  process  involving  relatively  slow  chemical  reactions. 

In  order  to  utilize  the  observations  of  soot  production  in  a  stagnation  point  flame  and 
previous  analyses  of  reacting  stagnation  point  flow  s,  the  configuration  of  a  laminar  counter-flow 
mixing  layer  has  been  adopted  for  the  present  numerical  study.  Specifically,  the  flow  under 
investigation  is  in  the  vicinity  of  the  forward  stagnation  point  of  a  porous  cylinder  from  which 
a  gaseous  reactant  is  ejected.  The  other  reactant  is  in  the  flow  which  approaches  the  cylinder. 

The  balance  of  reaction  and  convectivc-dilfusive  transport  in  this  flow  is  alfceted  by  the  rate 
of  strain,  i.e.  the  velocity  gradient.  The  phcnonienon  of  flame  stretch  has  been  shown  to  result 
in  the  extinction  of  reaction  in  a  stagnation  point  dilliision  flame  (Fcndell,  1965,  Linan,  1974). 
In  the  present  study  the  effect  of  the  strain  rale  t>r  vchK'ity  gradient  on  nucleation  kinetics  is 
investigated  for  two  generalized  isothermal  cases,  one  in  which  the  formation  of  condensable 
monomer  is  rate  limiting,  the  other  for  which  the  subsequent  growth  reactions  arc  rate  limiting. 
It  will  be  shown  for  both  cases  that  when  the  time  scale  for  the  rate  limiting  process  is  comparable 
to  the  time  scale  for  mixing  by  convection  and  diflusion,  the  fluid  mechanics  can  have  a  significant 
effect  on  the  nucleation  process. 


2.  HOMOGENEOUS  NUCLEATION 


2.1  Steady  State  Theories 

A  saturated  vapour  does  not  condense  to  the  liquid  phase  in  the  absence  of  foreign  nuclei 
until  a  critical  supersaturation  is  reached.  At  this  condition  the  onset  of  nucleation  occurs  and 
the  supersaturated  vapour  quickly  condenses  into  liquid.  Condensation  of  the  supersaturated 
vapour  occurs  by  the  formation  of  clusters  of  monomer  units  which  are  drawn  together  by 
intermolecular  attractive  forces.  The  build-up  of  these  clusters  can  be  represented  by  a  process 
in  which 


Cl  -t-  Ci_i  -.-Cj  (./=2, . M) 


(1) 


The  increase  in  size  of  a  cluster  Cj  i  is  due  to  the  incorporation  of  a  single  monomer.  Collisions 
of  clusters  with  other  clusters  are  assumed  to  be  relatively  unimportant  in  classical  theories 
because  of  the  predominance  of  monomer  in  the  mixture.  Recent  theoretical  evidence  suggests 
that  this  assumption  might  not  always  be  justified  (Gelbard  and  Seinfeld,  1979)  for  low  super¬ 
saturations.  Clusters  can  also  be  reduced  in  size  by  the  evaporation  of  monomer  in  the  reverse 
process. 

In  order  to  describe  the  concentration  of  clusters  of  various  sizes  in  a  system  the  classical 
theory  resorts  to  identification  of  cluster  properties  with  those  of  the  condensed  bulk  phase. 
This  is  the  so-called  capillarity  assumption.  Clusters  are  treated  as  small  liquid  droplets  with  a 
bulk  free  energy  and  a  surface  free  energy.  Formation  ofa  cluster  from  the  vapour  phase  monomer 
involves  a  decrease  in  bulk  free  energy  for  the  association  of  the  monomer  but  an  increase  in 
the  surface  free  energy.  The  change  in  the  Gibbs  free  energy  can  be  written  in  terms  of  liquid 
phase  properties 
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in  which  S  is  the  supersaturation  ratio  (T'/f’.-ii)  and  a  is  the  surface  tension  for  a  flat  liquid  surface. 
If  ACn  is  plotted  as  a  function  of  the  number  of  monomer  units  in  a  cluster  for  a  supersaturated 
vapour  then  a  maximum  is  evident  at  n*  (figure  I).  For  this  cluster  the  increase  in  AC„  due  to 
surface  free  energy  is  counterbalanced  by  the  decrease  in  bulk  free  energy  due  to  association 
of  n*  monomers.  A  cluster  with  n*  monomer  units  is  in  a  metastable  equilibrium  with  the  vapour; 
evaporation  occurs  at  the  same  rate  as  accretion  of  monomer.  This  cluster  is  referred  to  as  the 
critical  size  cluster  or  nucleus  and  the  rate  of  formation  of  clusters  of  the  critical  izc  represents 
the  nucleation  rate.  Once  this  thermodynamic  barrier  is  surmounted  then  very  rapid  growth 
of  a  cluster  is  possible  with  the  subsequent  formation  of  macroscopic  condensate  particles. 

As  an  example  of  the  classical  theory,  Frenkel's  steady  stale  rate  for  nucleation  is  expressed 
in  terms  of  bulk-phase  properties  (Frenkel.  1955), 
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(3) 


In  this  equation  ^  is  a  mass  accommodation  or  sticking  factor  and  represents  the  fraction 
of  collisions  which  are  successful  in  adding  a  monomer  to  a  growing  cluster.  The  first  term  in 
brackets  represents  the  rate  of  accretion  of  monomer.  The  second  term  in  brackets  represents  the 
Zeldovich  non-equilibrium  factor  which  accounts  for  the  difl'ercnce  between  a  metastable  equili¬ 
brium  size  distribution  of  sub-critical  clusters  and  the  size  distribution  attained  under  steady 
state  conditions.  The  sensitivity  of  the  nucleation  rate  to  supersaturation  arises  through  the 
dependence  of  the  critical  cluster  radius  (r*)  on  supersaturation  (.5). 

/•*  iTYfJkThtS.  (4) 

Modern  theories  of  nucleation  have  extended  the  classical  theory  by  including  additional 
terms  in  the  expression  for  the  change  in  Gibbs  free  energy  with  cluster  formation.  Evaluation 
of  rotational  and  translational  terms  from  statistical  mechanics  is  relatively  straightforward 
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but  considerable  controversy  has  been  generated  over  the  so-called  replacement  free  energy 
(Zettlemoyer,  1969). 


2.2  Transient  Nucleation 

In  a  uniform  vapour  which  is  brought  to  a  state  of  critical  supersaturation  there  is  a  time 
delay  until  the  onset  of  steady  state  nucleation.  For  most  simple  substances  the  estimated  time 
lag  is  of  the  order  of  1  /x  sec  which  is  insignificant  in  terms  of  time  scales  for  many  practical 
situations.  When  chemical  reactions  are  involved  in  nucleation,  however,  the  time  lag  may  be 
lengthened  so  that  transient  effects  are  significant  throughout  a  substantial  part  of  a  system  as 
suggested  by  Courtney  (1963)  for  a  reacting,  condensing  rocket  exhaust. 

Numerical  calculations  (Courtney  (1962),  Abraham  (1969),  Bauer  and  Frurip  (1977))  of 
the  approach  to  steady  state  involve  integration  of  a  master  equation  for  the  population  of  clusters 
in  the  system.  The  transition  probabilities  can  be  based  on  a  droplet  model  with  the  use  of  bulk- 
phase  properties  as  in  equation  (2).  Alternatively,  a  Gibbs  function  defined  by  extensive  thermo¬ 
dynamic  quantities  can  be  used  to  relate  forward  to  backward  transition  rates,  requiring  the 
evaluation  of  the  change  in  standard  enthalpy  (A//^)  and  entropy  (AA„)  for  the  association  of 
n  monomer  units  into  a  cluster.  This  is  the  approach  developed  by  Frurip  and  Bauer  (1977). 
It  involves  the  use  of  a  semi-empirical  formulation  of  the  enthalpy  and  entropy  changes  but 
avoids  the  use  of  the  capillarity  assumption  which  is  highly  dubious  for  small  clusters. 

Freund  and  Bauer  (1977)  examined  theoretical  and  experimental  results  for  the  energy  of 
condensation  of  an  atom  onto  a  cluster.  All  the  substances  were  metals  with  the  exception  of 
argon.  They  found  that  the  data  could  be  fitted  reasonably  well  by  t'lc  expression 

A£/// =  A£^  (1  -w  0--25)  (5) 

in  which  A£/«  is  the  energy  of  condensation  per  atom  for  it  vapour  atoms.  A£,  is  the  energy  of 
vapourization  or  sublimation  of  the  bulk  phase,  and  «  is  the  number  of  atoms  in  the  cluster. 
The  binding  energy  per  atom  in  a  cluster  is  less  than  in  the  bulk  condensed  phase  as  reflected 
in  the  lower  average  density  of  a  cluster.  In  order  to  calculate  the  Gibbs  runction  the  enthalpy 
change  for  the  association  of  it  atoms  or  molecules  is  assumed  to  be  approximately  equal  to 
the  energy  change,  i.e. 

A//n  ~  A£n. 

Under  the  conditions  of  interest  this  approximation  will  introduce  errors  of  the  order  of  10-20";, 
in  A/7„. 

Calculation  of  the  Gibbs  function  for  the  association  of  ii  monomers  is  completed  by  the 
estimation  of  the  entropy  change  involved  in  this  process.  Bauer  and  Frurip  (1977)  analysed 
the  entropy  change  in  terms  of  changes  in  a  structural  entropy  component  (incorporating  trans¬ 
lational,  rotational  and  vibrational  contributions)  and  a  configurational  term  which  accounts 
for  the  number  of  iso-energctic  arrangements  that  can  be  adopted  by  a  cluster  of  monomers. 
Consideration  of  experimental  data  for  polymerization  reactions  and  the  limiting  entropy  change 
for  bulk  condensation  lead  Bauer  and  Frurip  to  an  empirical  expression  for  the  structural  entropy 
decrease  in  the  association  of  it  monomers,  viz.; 

AS,,//;  -29n-^/(3  +  i/-)  (6) 

The  configurational  term  is  estimated  by  counting  the  number  of  isomers  which  arc  possible  for 
small  clusters  together  with  a  matching  algebraic  expression  for  larger  clusters. 

The  standard  Gibbs  free  energy  change  for  the  association  of  it  vapour  phase  monomers 
into  a  cluster,  viz.:  it/(  ->  An  is 

AC„  =  A//„  7-AS„.  (7) 

For  the  addition  of  monomer  to  a  cluster  of  size  ii.  the  free  energy  increment  for  the  process 


yields  the  connection  between  the  forward  and  backward  rate  constants 


(AO„,i  AGn) 
RT 


(8) 


in  which  A  has  units  of  [m'*  atm  mole  *  K  *1-  With  the  ratio  between  Au  and  A.  u  fixed  by  equation 
(8)  it  is  then  possible  to  integrate  the  master  equation  for  the  population  of  different  size  clusters 
in  the  system. 

Bauer  and  Frurip  (1977)  studied  the  Fe  vapour  system  at  1600  K.  They  found  that  steady- 
state  conditions  were  reached  in  the  order  of  10  '  secs.  Their  approach  yielded  a  criterion  for  a 
cluster  of  critical  size  in  kinetic  terms  as  one  for  which  the  number  flu.xes  of  clusters  into  and  out 
of  that  size  group  were  approximately  equal.  The  presence  of  a  kinetically  determined  critical 
cluster  was  found  not  to  require  a  maximum  in  the  Gibbs  function  as  it  does  in  the  classical 
theory. 

The  kinetic  approach  does  not  rely  on  the  use  of  an  equilibrium  size  distribution  w  ith  a  steady 
state  correction  as  in  the  classical  analysis.  Unfortunately  the  method  is  purely  computational 
and  fails  to  provide  an  analytical  solution  to  the  problem.  At  present,  however,  it  seems  that  a 
numerical  approach  is  necessary  if  one  is  to  consider  the  entire  history  of  a  condensing  vapour 
with  the  inclusion  of  the  initial  transient  phase. 


2.3  Chemical  Nucleation 

Nucleation  of  condensable  ^  ape.  .  s  in  a  number  of  reacting  flows  can  exhibit  transient-like 
effects.  Courtney  (1963)  conchf-j.  i  that  the  condensation  of  some  non-volatile  reaction  products 
in  rocket  exhausts  fueled  by  lithium  and  boron  propellants  could  be  described,  in  some  circum¬ 
stances,  as  a  non-steady  state  nucleation  process  because  the  time  lag  was  comparable  to  the 
residence  time  in  a  nozzle.  The  room  temperature  nucleation  of  NH  iCI  particles  from  reacting 
HCI  and  NHs  vapours  has  been  studied  by  Henry,  Gonzalez  and  Peters  (1983)  in  a  flow  reactor. 
They  reported  induction  times  of  the  order  of  2  seconds  between  the  initial  rapid  mixing  of  the 
two  reactants  and  the  onset  of  nucleation.  This  period  of  the  flow  is  certainly  non-steady.  These 
two  examples  illustrate  situations  in  which  the  production  of  condensable  monomer  is  slow  relative 
to  time  scales  of  a  flow.  The  production  of  soot  in  a  flame  is  a  case  in  w  hich  the  nucleation  process 
itself  is  relatively  slow.  The  effects  of  mixing  and  fluid  flow  on  nucleation  in  these  systems  has 
been  given  little  attention  in  the  past. 

A  generalised  model  has  been  set  up  to  investigate  nucleation  phenomena  in  a  counter-flow 
stagnation  point  configuration.  A  condensable  monomer  is  assumed  to  be  produced  by  the 
irreversible  one-step  reaction  between  reactants  A  and  B  wliich  are  diluted  to  initial  mass  fractions 
of  O'  1  in  an  inert  carrier  gas.  Reactant  A  is  in  the  gas  stream  issuing  from  a  porous  cylinder 
and  reactant  B  is  in  the  free-stream  flow  as  .shown  in  figure  2.  The  reaction  between  the  two 
species  is  assumed  to  follow  tlie  stoichiometry  described  in  equation  (9) 

zl  r  e  P  (9) 

in  which  the  combination  of  I  mole  of  A  and  i  mole  of  B  combine  to  yield  1  mole  of  product  P. 
This  stoichiometry  locates  the  zone  of  maximum  formation  rate  of  product  at  the  stagnation 
point  of  the  flow.  It  will  be  seen  later  that  this  choice  of  stoichiometry  facilitates  interpretation 
of  the  results  without  a  concomitant  reduction  in  generality. 

As  product  is  formed  it  is  allowed  to  simultaneously  nucleate  by  the  formation  of  clusters 
via  reversible  reactions  of  the  form 

C,  f  Gi  C, 

Gi  r  f  n  I  (  n  3, . A/. 

Two  limiting  cases  have  been  studied.  Tliey  correspond  to  firstly  the  situation  in  which 
production  of  condensable  monomer  is  slow  compared  willi  the  subsequent  nucleation  steps 
and.  secondly,  to  the  case  w  hen  the  nucleation  steps  arc  slow  compared  w  ith  monomer  production 
rates  because  of  small  accommodation  factors.  The  latter  situation  can  arise  when  cluster 
growth  involves  chemical  reaction  with  an  activation  energy  which  can  be  associated  with  a 
non-unity  accommodation  factor. 
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Both  cases  have  been  studied  for  isothermal  conditions.  Rates  of  monomer  production 
for  the  first  case  have  been  chosen  to  correspond  approximately  to  the  reaction  of  NH3  and  HCI 
at  400  K  with  a  rate  constant  of  approximately  10  to  100  m*/s-mole  (Countess  and  Heicklen, 
1973).  The  standard  Gibbs  free  energy  function  has  been  determined  from  the  heat  of  sublimation 
of  NHiCl  (39-9  kcahmole)  by  the  approach  of  Bauer  and  Frurip  ( 1977)  discussed  in  Section  2.2. 
Monomer-cluster  collision  rates  are  based  on  hard  sphere  rates  from  kinetic  theory  and  the 
accommodation  factor  for  monomer-cluster  collisions  is  unity.  Use  of  these  numerical  values 
is  only  a  convenient  means  for  ascribing  the  rate  constants  and  Gibbs  free  energy  functions  in 
this  purely  numerical  study;  the  use  of  different  Gibbs  functions  demonstrates  the  generality 
of  the  analysis.  The  rate  constants  occur  only  in  combination  with  the  flow  velocity  gradient 
(as  a  Damkohler  number  in  Section  4)  and  their  absolute  values  are  not  of  immediate  concern 
in  this  study.  Accurate  values  of  the  rate  constants  would  be  necessary  if  an  attempt  were  made 
to  deduce  velocity  gradients  from  a  Damkohler  number  for  a  specific  case. 

In  the  second  case  monomer  production  and  monomer-cluster  collisions  are  again  taken  to 
occur  at  hard  sphere  collision  rates  derived  from  kinetic  theory.  However,  in  this  case  the  growth 
rate  of  clusters  is  slow  compared  with  the  production  rate  of  condensable  monomer  due  to  a 
small  accommodation  factor.  This  factor  accounts  for  the  fraction  of  monomer-cluster  collisions 
which  are  successful  in  leading  to  incorporation  of  the  monomer  into  the  cluster.  For  many 
simple  substances  it  is  of  order  unity.  The  situation  of  present  interest,  however,  describes  a 
nucleation  process  in  which  the  accommodation  factor  is  much  less  than  one  due  to  the  involve¬ 
ment  of  chemical  reactions  in  the  bonding  of  monomer  to  cluster.  The  accommodation  factor 
in  the  calculations  has  been  chosen  as  10  '.  At  400  K  this  value  corresponds  to  an  activation 
energy  for  monomer  bonding  of  7-4  kcal,  mole  and  at  20(X)  K  to  an  activation  energy  of  36  kcal/ 
mole.  The  latter  value  is  repre.senlative  of  reaction  activation  energies  in  a  difl'usion  flame. 


3.  STAGNATION  POINT  FLOW 


Figure  2  illustrates  the  flow  configuration  that  has  been  studied.  Gas  which  is  seeded  with 
reactant  A  is  ejected  uniformly  around  the  circumference  of  a  porous  cylinder.  The  cylinder 
is  situated  in  an  approaching  flow  which  is  seeded  with  reactant  and  which  is  of  tmiform 
velocity  far  upstream  of  the  cylinder.  The  region  of  interest  is  in  the  vicinity  of  the  forward 
stagnation  point  where  a  local  similarity  solution  for  the  flow  field  is  valid.  The  advantage  of 
analysing  this  stagnation  point  flow  is  that  a  similarity  solution  permits  the  reduction  of  the 
partial  dilTerential  equations  for  mass,  moincntum  and  energy  to  ordinary  difl'erential  equations 
with  one  independent  variable. 

The  boundary  layer  form  of  the  momentum  equation  in  the  .v-direction  is  given  by 


(Vn  a;/  a  I  .'u\  (//’ 

'"'s  xU's  / 

a\  \  I  <l\ 


In  the  invi.scid  flow  outside  the  bouiulary  layer  the  pressure  gradient  is 


(//’ 

^/v 


/'  " 


t/n,- 

r/.v 


(10) 


(II) 


P  II  a 

where  a  is  the  imposed  velocity  gradient. 

A  transformed  normal  coordinate  is  introduced 


O' 


(12) 


The  equation  of  continuity  is  satisfied  by  relating  both  velocity  components  to  a  modified  stream 
function,  /,  so  that 


U  «r  /  ' 

where  the  prime  denotes  a  derivative  with  respect  to  >/  and 

i-  (m  )'  -  / 


(13) 


(14) 


6 


After  some  manipulation  the  transformed  momentum  equation  becomes 
/"' +#"-r^  + 1  =  0.  (15) 

The  boundary  layer  for  the  conservation  of  mass  fraction,  Tj,  of  species,  /,  is 
3  Ki  D  /  ^  ^ 

pu-  +pv  =  -  -WiWi.  (16) 

In  this  equation  is  a  binary  diffusivity  of  species  i  in  the  mixture.  The  reaction  rate  term,  Wi, 
is  the  net  molar  rate  of  production  of  species  /  where  nn  is  the  molecular  weight  of  this  species 
and  the  stoichiometric  coefficient  for  the  reaction  of  equation  (9)  is  unity. 

When  the  similarity  transformation  is  applied  to  this  equation  the  result  is 


1 

Sci 


y"+fYi' 


WiVVi 

ap 


where  (Schmidt  number)  for  species  /  is  defined  as 


(17) 


For  an  irreversible  one-step  reaction  such  as 

A  +  B  ~>P 

the  rate  of  production  of  product  P  is  expressed  as 
wp  =  kv  [A]  [S] 

In  terms  of  mass  fractions  equation  (20)  reads  as 

,  pY\  pYft 
Wp  —  kp 

ntA  niK 

The  conservation  equation  for  the  mass  fraction  of  product  P  is 


'  Yp"+fYp'  =  (  Ka  Yl^. 

Scp  \m,\  ttitt  a  ) 


The  first  group  of  terms  on  the  right  hand  size  defines  a  Damkohler  number 

/HP  pk 

Uv  =  — 

m\  m\\  a 


(18) 

(19) 

(20) 

(21) 

(22) 

(23) 


which  expresses  a  ratio  between  a  characteristic  flow  time  (through  a  ^  and  a  characteristic 

chemical  reaction  time  ( through  *  Y  Boundary  conditions  at  the  cylinder  surface  and 

V  mp  pk ) 

in  the  free  stream  are  required  to  complete  the  specification  of  the  problem. 

At  the  cylinder  wall  the  ejection  velocity  of  gas  is  cw.  The  stream-function  at  the  wall 
(tj  =  0)  is  then  given  by 

/(O)  =  —  rw/ v'o*'-  (24) 

Reactant  A  is  ejected  from  the  cylinder  and  reactant  B  is  in  the  frec-stream  flow  (see  figure  2). 
A  convective-diffusive  balance  at  the  wall  gives  the  follosving  boundary  conditions  for  the  mass 
fractions  of  A  and  B 

Ta(0)  -  Yi  -  Fa'(0)/S(  a  ./(O)  (25) 

and 

Fn(0)-  Tn'(0)/5:c„./(0),  (26) 

Yi  is  the  mass  fraction  of  inert  carrier  gas. 
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In  the  free-stream 


./■  =  I. 

=  0, 

Kii  =.  I  -  Yi. 


and 


(27) 


Nucleation  is  treated  by  solving  a  set  of  simultaneous  equations  of  the  form  of  equation  (22) 
for  the  mass  fractions  of  clusters.  The  Schmidt  number  of  a  cluster  is  derived  from  the  kinetic 
theory  for  diffusivities  (Hirschfelder,  Curtiss  and  Bird  (1954)).  The  rate  constant  for  a  collision 
between  a  cluster  with  /  monomer  units  and  a  single  monomer  is  assumed  to  increase  with  as 
spherical  cluster  sizes  increase. 

Solution  of  the  problem  entails  the  definition  of  the  flow  field  through  equation  (15)  for  a 
given  ejection  velocity  and  approach  velocity  followed  by  solution  of  the  coupled  equations  for 
Ka  and  Kb-  This  step  is  followed  by  the  simultaneous  solution  of  coupled  equations  for  the  pro¬ 
duction  of  monomer  and  larger  clusters  with  kinetics  (discussed  in  Section  2.3)  that  are  character¬ 
ised  by  a  Damkohler  number  for  monomer  production  (£>i)  and  Damkohler  numbers  (Di, 

i  =  2 . .  W)  based  on  the  rate  constant  and  molecular  weight  for  the  growth  of  a  cluster 

with  /■  monomer  units.  The  results  are  reported  for  conditions  specified  by  D\  and  Dz  where 
Di  is  the  Damkohler  number  for  monomer-monomer  collisions. 


4.  SOLUTION  PROCEDURE 


The  equation  for  momentum  (15)  and  the  equations  for  species  mass  fractions  (22)  are  written 
in  central  difference  forms.  The  solutions  of  the  resulting  tri-diagonal  systems  of  equations  are 
obtained  implicitly  with  a  Gauss  elimination  procedure  for  120  equally  spaced  grid-points  with 
an  interval  in  ij  of  0  04.  Iteration  of  the  calculations  is  continued  until  the  relative  difference 
between  successive  iterations  at  a  grid  point  is  less  than  0  001.  The  non-linearities  in  the 
momentum  equation  and  the  equation  for  monomer  mass  fraction  (the  term  involving  Ki  x  Ki) 
are  handled  by  quasi-linearization  (Bellman  and  Kalaba,  1965).  Refinement  of  the  grid  spacing 
by  25%  to  intervals  in  i/  of  0  03  results  in  differences  in  computed  cluster  concentrations  of  less 
than  1  %. 

The  non-dimensional  blowing  rate  from  the  cylinder  is  /(O)  —  10  for  all  cases.  The 

stoichiometry  of  reaction  (equation  (9))  between  reactants  in  the  flow  from  the  cylinder  and  in 
the  free-stream  is  chosen  so  that  the  stoichiometric  coefficients  of  A  and  B  in  (9)  are  I,  i.e.  one 
mole  of  A  reacts  completely  with  I  mole  of  B.  Both  reactants  are  diluted  in  an  inert  carrier  gas 
so  that  their  mass  fractions  in  the  feed-streams  are  OOOl  and  the  mass  fraction  of  inert  gas  is 
constant  throughout  the  flow. 

Damkohler  numbers  are  based  on  a  gas  density  of  I  •  18  kg  m  •*.  stoichiometric  coefficients 
of  1,  and  a  molecular  weight  of  condensable  monomer  (NH  jCI)  of  53  kg  kmole  '.  A  Damkohler 
number  of  10*  with  a  bimolecular  collision  rate  constant  of  10*  m*  mole  ’  sec  '  corresponds 
to  a  velocity  gradient  of  23  -7  sec  ’.  If  the  cylinder  radius  were  10  mm  then  the  approach  flow 
velocity  which  is  required  far  upstream  to  achieve  this  condition  would  be  118  mm  sec  '.  It 
follows  that  the  approach  flow  velocity  which  is  required  to  yield  a  Damkohler  number  of  10® 
is  1 1  -  8  m  sec  *. 

Damkohler  numbers  for  the  production  of  clusters  arc  modified  by  increases  in  the  molecular 
weight  of  clusters  and  by  increases  in  the  collision  rate  as  cluster  size  increases.  For  clusters  of 
particles  which  incorporate  more  than  20  monomer  units  averaged  conversion  rates  arc  used 
in  the  same  manner  as  Jensen  (1974).  so  that  if  the  number  of  monomers  in  a  cluster  is  «i  and 
«i  <  H|  then  for  the  process  Ci  +  Ci  Cj 


</lQ] 

(It  dt 
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(28) 


The  Damkohler  number  D\  is  based  on  the  averaged  rate  constant  k\i(n\  n,). 

A  limit  must  be  set  in  the  calculations  on  the  maximum  allowable  size  of  cluster  or  particle. 
The  present  calculations  involve  removal  of  particles  greater  than  this  size  from  the  system. 
Because  cluster  cluster  interactions  are  not  considered,  the  primary  effect  of  changing  the  upper 
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size  limit  is  to  change  the  rate  of  loss  of  monomer  via  particle  growth.  All  the  results  which  are 
presented  are  obtained  by  considering  particles  containing  up  to  240  monomer  units.'  Some 
cases  have  been  computed  with  a  maximum  allowed  particle  size  of  10  monomer  units;  the 
effect  of  cluster  concentrations  compared  with  the  full  calculation  is  less  than  0  01%.  This 
is  due  to  very  low  concentration  of  clusters  which  contain  more  than  a  few  monomer  units; 
extension  of  the  computations  to  clusters  of  240  units  serves  essentially  to  give  a  clearer  represen¬ 
tation  of  the  behaviour  of  the  system.  However,  termination  of  the  calculations  at  240  monomers 
ignores  the  role  of  larger  aerosol  particles  which  can  contribute  to  significant  heterogeneous 
nucleation.  This  question  is  addressed  in  a  later  section. 


5.  EFFECT  OF  MIXING  ON  NUCLEATION 
5,1  Reacting  Stagnation  Point  Flow  Field 

The  reactions  which  arc  included  in  the  model  are  not  significantly  endothermic  or  exo¬ 
thermic  at  low  concentrations.  Consequently,  there  is  no  effect  of  chemical  reaction  or  conden¬ 
sation  on  the  flow  field.  The  flow  is  the  same  for  all  the  cases  which  have  been  modelled  and 
can  be  solved  independently  of  the  nucleation  kinetics.  At  the  cylinder  wall  the  non-dimensional 
blowing  velocity  /(O)  is  —1  0.  The  two  velocity  components,  viz.:  i’  =  and  the  tan¬ 

gential  component  ujue  where  We  is  the  velocity  at  the  boundary  layer  edge,  are  shown  in  figure  3 
as  functions  of  the  non-dimensional  normal  coordinate,  ri  where 


The  tangential  velocity,  a,  exhibits  a  laminar  boundary  layer  profile.  The  thickness  of  the 
layer  scales  with  (v/a)'^^.  The  normal  velocity  component  indicates  that  the  stagnation  point  is 
at  about  rj  =  1-9.  The  distance  of  the  stagnation  point  from  the  cylinder  will  vary  with 
(vjaY'^  for  a  fixed  gas  ejection  velocity.  Calculations  of  nucleation  kinetics  are  performed  with 
these  flow  conditions. 

The  interaction  between  the  flow  field  and  nucleation  is  determined  by  the  Damkohler 
number  as  discussed  in  Section  4.  If  the  Damkohler  number  is  small  then  chemical  reaction  is 
slow  compared  to  residence  time  in  the  flow.  The  effect  of  Damkohler  number  on  the  reaction 
of  two  species  contained  in  the  two  mixing  streams  is  shown  in  figure  4.  For  the  reaction  rates, 
densities  etc.,  cho.sen  in  Section  5,  a  Damkohler  number  of  10*  corresponds  to  a  velocity  gradient 
of  23  -7  sec"'.  With  one  mole  of  reactant  A  consuming  one  mole  of  reactant  B  the  reaction  zone 
is  located  in  the  vicinity  of  the  stagnation  point.  For  low  velocity  gradients  (£>i  =  10*)  the  region 
of  reaction  is  narrow  with  little  overlap  of  reactants.  As  the  velocity  gradient  is  increased  (Di  = 
10*)  there  is  considerably  more  interpenetration  of  reactants. 

5.2  Free  Energy  Functions 

Standard  Gibbs  free  energy  functions  for  use  in  the  computation  of  nucleation  kinetics 
have  been  determined  using  the  empirical  approach  of  Bauer  and  Frurip  (1977).  Enthalpies 
and  entropies  are  estimated  from  curves  fitted  to  theoretical  and  experimental  data  (see  Section 
2.2).  The  substances  considered  by  Bauer  and  Frurip  were  monatomic  species;  the  use  of  their 
empirical  correlations  for  polyatomic  species  involves  some  uncertainty.  However,  it  will  be 
shown  that  the  choice  of  the  Gibbs  function  docs  not  significantly  alter  the  overall  behaviour 
of  the  nucleation  kinetics.  The  system  that  has  been  modelled  in  most  cases  consists  of  NH.iCl 
vapour  for  which  the  estimated  Gibbs  functions  at  400  K  and  550  K  are  shown  in  figure  5  in  a 
dimensionless  form. 

At  400  K  the  Gibbs  function  very  quickly  attains  large  negative  values  which  indicate 
stable  clusters.  There  is  no  maximum  in  the  curve.  At  550  K  (61  K  below  the  temperature 
of  sublimation)  the  Gibbs  function  exhibits  a  maximum  at  a  cluster  with  3  molecules  and  does 

'  Based  on  the  bulk  density  of  NHiCI  the  diameter  of  these  particles  is  about  3  nm. 
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not  attain  negative  values  until  clusters  contain  more  than  6  molecules.  Also  shown  in  figure  5 
is  the  Gibbs  function  for  MgO  at  1700  K.  Neville  and  Sarofim  (1982)  reported  a  study  of  particle 
production  from  inorganic  vapours  evolved  during  the  combustion  of  coal.  Magnesium  oxide 
particles  were  formed  by  the  nucleation  of  the  oxide  vapour  which  was,  in  turn,  produced  by 
the  oxidation  of  magnesium  vapour  around  a  coal  particle.  The  magnesium  oxide  system  affords 
an  additional  case  to  which  a  kinetic  scheme  can  be  applied  and  permits  an  evaluation  of  the 
generality  of  the  conclusions  which  have  been  reached.  It  can  be  seen  from  figure  5  that  MgO 
vapour  condenses  quite  readily  at  1700K  with  the  Gibbs  function  rapidly  attaining  negative 
values.  Clusters  with  more  than  4  or  5  molecules  are  quite  stable  according  to  these  estimates 
of  the  free  energy. 


5.3  Nucleation 

5.3.1  Fast  Vapour  Formation:  Slow  Nucleation  (Di  ^  1,  C  10  “*) 

If  the  formation  of  a  cluster  involves  chemical  bonding  then  the  accommodation  or  sticking 
factor  for  monomer-cluster  interactions  will  be  less  than  one.  The  accommodation  factor,  is 
chosen  to  be  10^^  in  these  calculations  (see  Section  2.3).  Damkohler  numbers  for  cluster  forma¬ 
tion  are  then  related  to  D\  as  Z)j  =  {  J  =  2, . . M  where  the  terms  in  j  account 

for  variations  in  molecular  weights  and  collision  rates  with  cluster  size. 

In  the  absence  of  a  simple,  chemically  reacting  system  to  which  this  model  can  be  applied, 
the  Gibbs  function  for  NH  iCI  vapour  at  400  K  has  been  used  in  the  calculations.  This  choice 
aids  in  the  comparison  of  the  results  with  those  obtained  in  the  following  section. 

Three  velocity  gradient  conditions  have  been  used  to  yield  Di  of  10®,  10’  and  10®.  Calcu¬ 
lated  mass  fractions  for  various  size  clusters  or  particles  are  shown  in  figure  6  at  these  three  values 
of  Damkohler  number.  It  is  apparent  from  figure  6  that  the  nucleation  process  is  very  sensitive 
to  flow  conditions  when  cluster  growth  is  slow.  For  D\  =  10®,  the  mass  fraction  of  clusters  with 
80  molecules  is  2x  10-®  which  corresponds  to  about  3x  10*®  clusters  m“®.  When  the  approach 
flow  velocity  is  increased  by  a  factor  of  10  so  that  D\  =-  10’,  the  mass  fraction  of  clusters  with 
80  molecules  reduces  by  three  orders  of  magnitude  to  about  2x  10  '*  or  3x  10*'^  clusters  m  ®. 
With  a  further  increase  in  the  approach  flow  velocity  by  an  order  of  magnitude  so  that  Di  =  10®, 
the  concentrations  of  clusters  of  any  significant  size  are  negligible.  The  production  of  larger 
particles  has  been  suppressed,  in  effect,  by  controlling  the  flow  field. 

The  sensitive  response  of  particle  formation  to  changes  in  the  velocity  gradient  is  illustrated 
more  directly  in  figure  7  for  a  range  of  Damkohler  numbers.  The  mass  fractions  and  number 
concentrations  of  particles  with  200  monomer  units  at  the  stagnation  point  of  the  flow  are  shown 
as  functions  of  a  Damkohler  number.  For  large  Damkohler  numbers  there  is  a  relatively  weak 
dependence  of  particle  concentration  on  velocity  gradient.  However,  at  a  Damkohler  number  of 
around  5  Ox  10’  kinetic  effects  start  to  become  dominant  and  the  concentration  of  particles 
begins  to  decrease  rapidly.  At  a  Damkohler  number  of  about  2x  10®  there  are  effectively  no 
large  particles  in  the  flow. 

It  was  pointed  out  in  Section  5  that  these  calculations  involve  the  removal  of  particles 
with  more  than  240  monomer  units  from  the  system.  In  fact,  with  a  large  Damkohler  number 
there  will  be  a  substantial  amount  of  larger  aerosol  particles  present  in  the  flow  as  the  results 
presented  in  figure  6  suggest.  For  D\  of  10*  there  is  little  diminution  of  mass  fractions  of  clusters 
with  up  to  240  monomer  units. 

Friedlander  (1982)  demonstrated  in  a  numerical  calculation  that  an  aerosol  can  scavenge 
monomers  and  clusters  and  can  reduce  the  nucleation  rate  significantly.  Condensation  onto 
an  existing  particle  provides  an  alternative  route  in  contrast  to  homogeneous  nucleation  for  gas 
to  particle  conversion.  The  present  calculations  have  been  limited  to  240  monomer  units  by 
the  restriction  of  available  computing  services  and  by  the  numerical  difficulty  of  incorporating 
large  particles  into  the  study.  The  latter  difficulty  arises  as  the  diffusivitics  of  large  particles 
become  quite  small  and  the  Schmidt  numbers  become  correspondingly  large.  The  inverse 
Schmidt  number  multiplies  the  highest  order  term  in  equation  (17).  As  the  Schmidt  number 
becomes  large  the  problem  tends  to  become  singular  as  the  governing  equation  is  reduced  from 
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second  to  first  order.  The  numerical  difficulties  which  arise  from  considering  large  aerosol 
particles  in  the  mixing  layer  have  precluded  their  consideration  in  this  study. 

The  behaviour  of  the  nucleating  mixing  layer  will  be  unaffected  by  an  aerosol  at  high 
velocity  gradients  (small  Damkohler  numbers)  for  which  the  mass  fractions  or  large  clusters  are 
negligible.  For  low  velocity  gradients  the  nucleation  rate  will  be  reduced  and  the  cluster  size 
distribution  will  be  modified  by  the  presence  of  an  aerosol.  However,  it  is  expected  that  the 
qualitative  behaviour  of  the  system  will  be  similar. 

A  comparison  of  the  results  presented  in  figure  6  with  those  of  Bauer  and  Frurip  (1977) 
for  a  kinetic  nucleation  model  of  iron  vapour  condensation  reveals  similarities  between  the  two 
sets  of  calculations.  The  sensitive  response  of  nucleation  kinetics  in  the  stagnation  point  flow  to 
changes  in  the  velocity  field  is  analogous  to  the  transient  conditions  computed  by  Bauer  and 
Frurip  for  their  nucleating,  uniform,  iron  vapour  system.  The  population  of  large  clusters 
was  very  small  in  their  calculations  for  periods  of  10“i**  to  10"®  secs  after  the  initiation  of  cluster 
growth  at  a  supersaturation  ratio,  S,  of  3040.  By  10"’  secs  their  system  had  reached  a  steady- 
state  with  significant  levels  of  large  clusters.  The  behaviour  of  the  nucleating  stagnation  point 
flow  appears  to  be  directly  analogous  to  this  transient  nucleation  phenomenon. 

The  distribution  of  clusters  of  different  sizes  throughout  the  flow  is  shown  in  figure  8  for  a 
Damkohler  number  of  10®.  The  peaks  in  cluster  concentration  occur  at  the  same  location  close 
to  the  stagnation  point  of  the  flow.  At  this  point  the  effects  of  diffusion  and  convection  of  clusters 
are  small  because  of  the  near  zero  concentration  gradients  and  velocity.  As  a  result  it  is  possible 
to  estimate  the  mass  flux  from  one  size  group  to  another  without  a  significant  complication  due 
to  transport  processes. 

Using  a  value  of  the  rate  constant  for  monomer-monomer  interactions  of  10®  m®  mole"’ 
sec"’  and  ^  =  10-‘>  it  is  possible  to  derive  the  mass  flux  in  cluster  size  space  (y)  in  one  direction 
for  a  given  size  cluster,  i.e.  yi-,i+i  is  the  mass  flux  (kg  m'®  s"’)  for  conversion  of  clusters  with  i 
monomers  to  i-fl  monomers  while  yi*.n-r  is  the  mass  flux  in  the  reverse  direction.  Table  1 
presents  these  mass  fluxes  for  two  cases,  viz.;  Di  =  10®  and  Di  =  10®. 

When  Di  =  10®,  i.e.  when  the  approach  flow  velocity  is  low,  the  forward  and  backward 
fluxes  are  balanced  for  small  clusters.  At  the  cluster  size  associated  with  the  ‘knee’  in  the  mass 
fraction  curve  of  figure  6,  the  forward  and  backward  fluxes  are  nearly  balanced,  and,  in  addition, 
the  two  forward  fluxes  into  and  out  of  this  size  group  are  approximately  equal  (within  an  order 
of  magnitude).  For  the  Di  ==  10®  case  the  ‘knee’  is  at  about  n  =  3  so  that  y2-,3  ~  ~ 

y3-.4  ~  y3,-4.  Bauer  and  Frurip  (1977)  identified  this  condition  as  the  specification  of  a  kinetically 
determined  critical  cluster  size.  The  critical  cluster  marks  the  point  of  a  switch-over  in  the  kinetics; 
the  population  of  small  clusters  is  established  by  balanced  forward  and  backward  conversion 
rates  while  the  population  of  larger  clusters  is  maintained  by  nearly  equal  forward  conversion 
rates  with  negligible  backward  rates  (for  example  yiz-is  ~  yi8-.i»  and  yi7,-i8  ~  yi8,-i9  ~  0). 
As  noted  by  Bauer  and  Frurip  (1977)  the  presence  of  a  critical  cluster  in  the  kinetics  calculation 
does  not  require  a  maximum  in  the  Gibbs  function  as  it  does  in  the  classical  theory.  The  com¬ 
puted  behaviour  of  the  system  modelled  with  NH4C)  vapour  at  400  K  for  which  there  is  no 
maximum  (figure  5)  is  essentially  the  same  as  at  550  K  (figure  9)  for  which  a  maximum  in  the 
Gibbs  function  is  evident. 

For  a  higher  approach  flow  velocity  ( Dy  =  10®)  Table  I  indicates  that  the  forward  and  reverse 
fluxes  for  clusters  with  more  than  one  monomer  unit  are  never  equal,  i.e.  the  reversible  processes 
involved  in  cluster  growth  are  not  equilibrated  even  for  small  clusters.  At  larger  sizes  the  pre¬ 
dominant  flux  is  the  forward  one  but  it  continues  to  decrease  with  increasing  cluster  size.  From 
these  results  it  is  evident  that  small  cluster  concentrations  do  not  have  a  chance  to  equilibrate 
within  the  time  scale  imposed  by  the  velocity  gradient.  As  a  result  larger  clusters  cannot  estab¬ 
lish  significant  concentrations.  This  behaviour  is  directly  analogous  to  a  transient  nucleation 
phenomenon. 

In  order  to  demonstrate  that  the  form  of  the  results  are  not  specific  to  a  particular  choice 
of  the  free  energy  function,  the  calculations  for  Dy  ~  10*  and  Dy  =  10®  have  been  repeated 
with  an  estimated  free  energy  function  for  NH4CI  at  550  K.  The  mass  fractions  for  various  size 
clusters  are  shown  in  figure  9.  The  bend  or  knee  in  the  curves  is  shifted  to  larger  clusters  (n~9) 
and  the  concentrations  are  substantially  lower.  If  the  cluster  for  which  the  four  fluxes  (yi_i_i, 
yui^i,  yi-i+i,  yi,-ni)  are  approximately  equal  (i.e.  at  the  ‘knee’  in  the  mass  fraction  curve) 
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is  taken  to  be  a  critical  size  cluster  or  nucleus  then  the  nucleation  rate  for  Di  =  10®  at  400  K 
is  of  the  order  of  I0^“  m  ®  s  '.  At  550  K  the  nucleation  rate  is  of  the  order  of  10’ m  ®  s  *. 

Despite  the  different  concentrations  of  large  clusters,  the  response  of  the  nucleation  process 
to  changes  in  Damkohler  number  is  much  the  same  at  400  K  and  550  K.  At  D\  ^  10*  there  is 
only  a  very  gradual  decay  in  the  number  concentration  of  large  particles  with  size  while  at 
Di  =  10®  there  are  effectively  no  large  particles.  The  critical  range  of  Damkohler  numbers 
which  determines  the  cessation  of  nucleation  evidently  does  not  depend  strongly  on  the  free 
energy  function.  It  must  be  noted,  however,  that  at  higher  temperatures  the  rate  constants  for 
cluster  growth  will  be  larger  and,  consequently,  the  velocity  and  velocity  gradient  must  be  corre¬ 
spondingly  larger  to  achieve  a  certain  Damkohler  number.  Whenever  the  nucleation  Damkohler 

numbers,  D,,  i  =  2 . M,  are  sufficiently  small  due  to  a  small  accommodation  factor  the 

flow  field  effect  reported  here  is  possible. 

Although  the  model  which  has  been  developed  here  for  nucleation  in  a  stagnation  point 
flow  is  isothermal  and  represents  a  limiting  case  for  which  condensable  monomer  production 
is  fast,  it  affords  some  insight  into  the  response  of  a  sooting  difl'usion  flame  to  variations  in  the 
local  flow  field.  The  formation  of  soot  occurs  at  high  temperatures  around  1800  K  on  the  fuel 
side  of  a  flame.  There  is  no  identifiable  single  species,  i.e.  monomer  which  leads  directly  to  soot 
particles  via  polymerization  reactions  although  the  polyacetylene  mechanism  of  Homann  and 
Wagner  (1967)  relies  upon  such  a  route  up  to  the  point  of  aromatic  ring  formation.  The  formation 
of  polycyclic  aromatic  compounds  (building  blocks  for  soot  particles)  from  straight  chain 
aliphatic  fuels  poses  a  mechanistic  problem  for  kinetic  schemes  of  soot  formation  (Cole.  Bittner, 
Longwell  and  Howard,  1984).  However,  once  aromatic  compounds  such  as  benzene  or  toluene 
have  been  formed  in  the  flame,  growth  of  more  stable,  larger  polyaromatic  species  can  proceed 
by  the  addition  of  abundant  stable  species  such  as  acetylene  or  highly  reactive  radicals.  Stein 
(1978)  has  used  group  additivity  methods  to  analyse  the  high  temperature  equilibria  of  acetylene 
and  polyaromatic  species.  He  suggested  a  route  which  is  favoured  thermodynamically  for  the 
formation  of  a  large  polycyclic  aromatic,  C<MiH.>|.  The  Gibbs  functions  for  the  formation  of 
these  polyaromatics  are  similar  in  form  to  the  Gibbs  functions  for  NH  iCI  association  which 
are  used  in  the  present  model. 

The  formation  of  condensable  monomer  is  assumed  to  be  much  faster  than  the  subsequent 
nucleation  in  this  model;  the  formation  of  acetylene  and  other  pyrolysis  products  in  a  flame  may 
not  approach  this  limiting  rate.  However,  data  for  stagnation  point  diffusion  flames  (Tsuji, 
1982)  indicate  constant  flame  temperatures  and  presumably,  constant  pyrolysis  rates  over  the 
range  of  velocity  gradients  in  which  soot  formation  is  suppressed. 

Despite  the  differences  between  the  model  and  the  complex  kinetics  of  soot  formation  in  a 
diffusion  flame  the  overall  nucleation  mechanisms  are  sufficiently  similar  to  allow  some  conclu¬ 
sions  to  be  drawn.  If  the  formation  of  polyaromatic  compounds  in  a  flame  involves  relatively 
slow  kinetics  at  some  stage,  a  large  velocity  gradient  in  a  stagnation  flow  may  reduce  species 
concentrations  substantially  below  equilibrium  levels.  The  result  of  non-equilibrium  concen¬ 
trations  of  the  least  stable,  small  preeursor  molecules  in  the  nucleation  chain  would  be  greatly 
reduced  rates  of  macroscopic  particle  formation. 


5.3.2  Slow  Vapour  Formation:  Fast  Nucleation  (l)i  ^  I,  Di  I,  i  2. . . .,  M) 

For  most  condensing  vapours  the  accommodation  cocflicicnt  t  ^  I  so  that  cluster  growth 
is  limited  by  the  collision  rate  between  monomer  and  clusters.  However,  when  a  condensable 
vapour  is  the  product  of  a  chemical  reaction  such  as  in  a  rocket  exhaust  (Courtney  (1963)), 
during  coal  combustion  (Katz  and  Donohue  (1982),  Neville  and  Sarolim  ( 1982))  or  a  simple 
room  temperature  reaction  such  as  HC  I  and  NH:i  (Countess  and  Heicklen  (1973),  Henry  el  al. 
(1983))  the  rate  of  monomer  formation  may  be  the  limiting  step  in  the  nucleation  of  particles. 
The  Damkohler  number  for  monomer  formation,  />i,  may  be  of  order  I  for  practical  reaction 
rates  and  velocity  gradients. 

Calculations  have  been  performed  for  the  HCI  NH:i  system  in  a  stagnation  point  mixing 
layer.  The  rate  constant  of  Countess  and  Heicklen  (1973)  for  the  HC  I  Nila  reaction  suggests 
values  of  D\  60.  6  and  ()-6  for  realistic  velocity  gradients.  1  he  computer  simulation  has  been 
carried  out  to  model  the  effect  of  varying  velocity  gradients  when  /)|  is  small.  As  l)[  decreases 


from  60  to  0-6,  Dz  for  monomer-monomer  collisions  decreases  from  10*  to  10®.  The  mass 
accommodation  factor  t  is  1.  The  Gibbs  free  energy  function  for  NH4CI  at  400  K  has  been  used. 

The  mass  fractions  of  clusters  at  the  stagnation  point  are  shown  in  figure  10  for  three  cases 
with  different  velocity  gradients.  Monomer  concentrations  decrease  by  two  orders  of  magnitude 
as  the  velocity  gradient  is  increased  by  two  orders  of  magnitude.  There  is  a  much  more  pronounced 
effect  on  the  concentration  of  larger  clusters.  An  increase  in  the  velocity  gradient  by  two  orders 
of  magnitude  results  in  13  orders  of  magnitude  decrease  in  large  cluster  mass  fractions.  The 
size  of  the  kinetically  determined  ‘critical’  cluster  increases  from  about  4  molecules  at  the  highest 
Damkohler  number  to  about  6  at  the  lowest  Damkohler  number.  Based  on  the  concentration 
of  these  clusters  the  nucleation  for  f)i  =  60  is  about  lO^*  m"*  s  '  and  for  Di  =  0-6  is  about 
10“*  m"*  s“i.  The  latter  rate  is  so  low  that  particle  formation  is  negligible. 

The  mass  fluxes  for  various  processes  are  shown  in  Table  2.  It  can  be  seen  that  forward 
and  backward  rates  are  equilibrated  for  small  clusters  for  both  high  and  low  Damkohler  numbers. 
This  behaviour  is  in  contrast  to  the  previous  case  (Section  5.3.1)  in  which  nucleation  was  slow 
compared  with  monomer  formation.  For  those  conditions  a  low  Damkohler  number  (high 
velocity  gradient)  prevented  equilibration  of  forward  and  baekward  rates  for  even  the  smallest 
clusters.  As  a  result  nucleation  behaviour  was  characteristic  of  a  transient  phenomenon.  In 
the  present  case  for  which  monomer  production  is  rate  limiting  the  behaviour  is  characteristic 
of  a  steady  state  system  for  all  the  velocity  conditions  which  have  been  computed.  The  response 
of  the  nucleation  process  to  velocity  gradient  changes  is  typical  of  the  sensitive  response  of  a 
nucleating,  uniform  vapour  to  variations  in  supersaturation. 

The  monomer  mass  fraction  when  Di  =  60  is  6x  10  ®  and  when  Di  ~  0  6  is  7x  10'’. 
Two  orders  of  magnitude  variation  in  supersaturation  can  give  rise  to  extremely  large  differences 
in  nucleation  rates  as  the  classical  theory  indicates  via  equations  (3)  and  (4).  Hence  the  effect 
of  flow  conditions  on  nucleation  rates  in  the  mixing  layer  for  this  case  can  be  interpreted  in  terms 
of  the  analogy  with  the  very  strong  effect  of  supersaturation  in  a  uniform  vapour.  The  influence 
of  flow  conditions  is  shown  directly  in  figure  1 1  for  the  mass  fractions  of  particles  with  200 
monomer  units.  The  curve  is  consistent  with  the  supersaturation  analogy. 

In  order  to  check  on  the  generality  of  these  observations  for  other  systems  a  calculation 
has  been  performed  with  the  estimated  Gibbs  function  for  MgO  vapour  at  1700  K.  The  Gibbs 
function  is  shown  in  figure  5.  It  corresponds  to  the  conditions  studied  by  Neville  and  Sarofim 
(1982)  for  MgO  particle  formation  during  coal  combustion.  In  keeping  with  the  findings  of  Neville 
and  Sarofim  (1982)  the  condensable  MgO  vapour  is  assumed  to  be  formed  by  the  gas-phase 
oxidation  of  magnesium  vapour.  Varying  the  Damkohler  numbers  by  changing  the  velocity 
gradient  over  two  orders  of  magnitude  produces  the  same  effect  for  the  MgO  system  (see  figure  12) 
as  the  NH4CI  system.  That  is,  mass  fractions  of  small  clusters  are  equilibrated;  the  response 
to  the  reduced  formation  rate  of  vapour  (reflected  in  two  orders  of  magnitude  decrease  in  monomer 
concentration)  is  equivalent  to  a  reduction  in  supersaturation  in  a  uniform  gas. 


CONCLUSIONS 

Chemical  nucleation  in  a  reacting  mixing  layer  has  been  modelled  numerically  with  an  iso¬ 
thermal,  stagnation  point  flow.  Reactants  in  two  mixing  streams  form  a  condensable  monomer 
product  which  subsequently  undergoes  nucleation  through  the  formation  of  clusters.  The  rele¬ 
vant  time  scale  for  convective-diffusive  transport  in  the  stagnation  point  flow  is  the  inverse 
rate-of-strain,  i.e.  the  inverse  velocity  gradient  in  the  vicinity  of  the  stagnation  point.  Calcula¬ 
tions  show  that  when  the  time  scale  for  the  rate  limiting  chemical  reaction  in  the  nucleation 
process  is  comparable  to  the  time  scale  imposed  by  the  flow  field,  nucleation  can  be  suppressed. 

Two  limiting  cases  have  been  investigated.  When  the  formation  rate  of  condensable  monomer 
is  slow,  nucleation  in  the  mixing  layer  exhibits  a  sensitive  response  to  the  velocity  gradient. 
The  concentrations  of  small  clusters  are  equilibrated  with  monomer  concentration  for  all  flow 
conditions  and  the  effect  of  the  flow  field  on  nucleation  is  analogous  with  the  effect  of  the  super¬ 
saturation  ratio  in  a  uniform  gas.  However,  when  monomer  cluster  reactions  are  rate  limiting 
because  of  a  small  accommodation  factor,  the  sensitive  response  of  the  nucleation  kinetics  to 
variations  in  velocity  gradient  follows  a  different  route.  Over  a  critical  range  of  velocity  gradients 
the  concentration  of  small  clusters  cannot  establish  equilibrium  levels.  The  non-equilibrium 
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concentrations  of  small  clusters  results  in  the  suppression  of  particle  formation  for  high  velocity 
gradients. 

The  suppression  of  soot  formation  in  a  stagnation  point  diffusion  flame  at  a  critical  velocity 
gradient  can  be  understood  in  terms  of  the  latter  mechanism.  The  effect  on  nucleation  kinetics 
of  the  mixing  rate  may  be  significant  in  other  reacting  flows  such  as  turbulent  diffusion  flames  and 
rocket  exhausts. 
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TABLE  1 


Mass  Fluxes  (kg  m~^s~0  Between  Cluster  Sizes 


Forward 


Di  =  108, 

c  =  10-4 

4-74x10-1 

4-74x10-1 

2-52x10-4 

2-26x10-4 

2-47x10-5 

1-18x10-5 

1-20x10-8 

1-77x10-9 

7-86x10-9 

3-3  X  10-81 

7-86x10-9 

3-3  X  10-81 

7-84x10-9 

3-3  X  10-81 

Reverse 


4-33x10-1 

4-24X  10-1 

2- 15  X  10-4 

9  -  58 X 10-5 

1-04x10-5 

2-51x10-’ 

Forward 


TABLE  2 

Mass  Fluxes  (kg  Between  Cluster  Sizes 

Di  =  60,  Z>2  =  10» 


1-32x10-8 

1  -32x  10-8 

4  -  77  X  10-13 

4-77x10-13 

3  -  55  X  10  15 

3-51  X  10-15 

2-43X  10  '« 

1  •99x  10-ie 

Di  =  0-6,  D 


Reverse 


1-76x10  18 

1  -  76  X  10-18 

7-32x  10  1" 

7-32X  10  19 

6  -  28  X  10-83 

6-28X  10-88 

5-02X  10  8« 

5  01  X  10  8« 

2  l2x  10  28 

2-08  X  10  88 
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STAGNATION  POINT  MIXING  LAYER  FLOW  USED  IN  THE  CALCULATIONS 
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JUTION  OF  CLUSTERS  THROUGH  THE  FLOW  FOR  CLUSTERS  WITH 
lONOMER  UNITS 


CLUSTER  MASS  FRACTIONS  AT  STAGNATION  POINT  FOR  SLOW  MONOMER 
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FIG.  12  CLUSTER  MASS  FRACTIONS  AT  STAGNATION  POINT  FOR  SLOW  MONOMER 
FORMATION  RATES  f  =  1  AND  AG^  FOR  MgO  AT  1700  K 
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the  nucleation  kinetics  is  analogous  to  the  transient  phase  of  nucleation  in  a  uniform  vapour. 


16.  Abstract  (Contd) 


17.  Imprint 

Aeronautical  Research  Laboratories,  Melbourne 


18.  Document  Series  and  Number  19.  Cost  Code  20.  Type  of  Report  and  Period  Covered 

Aero  Propulsion  Report  165  437410 


21.  Computer  Programs  Used 


22.  Establishment  File  Ref(s) 


